Exact Non-Markovian Master Equation and Dispersive Probing of Non-Markovian 

Process 
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For a bosonic (fermionic) open system in a bath with many bosons (fermions) modes, we derive 
the exact non-Markovian master equation in which the memory effect of the bath is reflected in the 
time dependent decay rates. In this approach, the reduced density operator is constructed from the 
formal solution of the corresponding Heisenberg equations. As an application of the exact master 
equation, we study the active probing of non-Markovianity of the quantum dissipation of a single 
boson mode of electromagnetic (EM) field in a cavity QED system. The non-Markovianity of the 
bath of the cavity is explicitly reflected by the atomic decoherence factor. 
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I. INTRODUCTION 



The open quantum system approach is of much sig- 
nificance due to its various applications in physics, e.g., 
quantum information, quantum transport, and quantum 
chemistry, etc. Since a realistic quantum system is in- 
evitably coupled to many degrees of freedom in its en- 
vironment that leads to decoherence of the systems, a 
general approach to the open quantum system is needed 
for its dissipative and dephasing process. The dynam- 
ics of an open system is conventionally described with 
three approaches: effective Hamiltonian quantum 
master equations 0, 0], and quantum Langevin equa- 
tions [H, The last two approaches are both based 
on the modeling with system plus bath, while the first 
one is phenomenologically given by a time-dependent or 
non-Hermitian Hamiltonian, which could lead to the dis- 
sipative motion equations. 

About twenty years ago, Yu and one (C. P. S.) of 
the authors revealed an intrinsic relation between the 
effective Hamiltonian and quantum Langevin equation, 
obtained from the Heisenberg equations [H, 0|. By 
discarding the quantum fluctuation for the wide wave 
packet, they derived the effective Hamiltonian of the sys- 
tem through the formally exact solution for the time- 
dependent wave function of the total system. However, 
the resulting effective Hamiltonian ignores the memory 
effect, which is induced by the back action of the bath 
with time delay. Therefore, if one wanted to recover 
the non-Markovian phenomenon with memory effect, the 
quantum fluctuation of the bath must be taken into ac- 
count in the above Heisenberg equation based approach. 
To this end, we need start from the Heisenberg equa- 
tions of the total system, which can reflect the original 
role of the bath. In this paper, without any approxima- 
tion, we derive the exact non-Markovian master equation 
of the system from the formal solution of the Heisenberg 



equations. The non-Markovian effect is contained in the 
time-dependent decay rates in a straightforward way (l4| . 

It is commonly believed that the Markov process hap- 
pens when the system-bath coupling is weak. However, 
with the rapid development of experimental technology, 
the strong-coupling limit can be reached. The theory of 
open quantum systems in the strong-coupling regime is 
required for a proper description of the non-Markovian 
dynamics. Recently, many works on exact quantum mas- 
ter equations have been done (Iol - [l8| . In particular, 
one (W. M. Z.) of the authors and his collaborators de- 
rived the exact non-Markovian master equations with a 
Lindblad-form for both Bose QJ, QJ] and Fermi Qj| 
systems by a path-integral method in coherent-state rep- 
resentation. We now revisit these non-Markovian mas- 
ter equations by generalizing our previous approach [f|, 
which was used to derive a partially factorized wave func- 
tion for open quantum system. Using the present gen- 
eralization to derive the reduced density matrix is quite 
straightforward. Here, we first construct the total den- 
sity matrix with the help of the formal solution of the 
Heisenberg equations, and then trace over the degrees of 
freedom of the bath to obtain the reduced density matrix 
of the system in the coherent state representation, in- 
stead of using the Feynman- Vernon influence functional, 
as was done in Refs. [l3l - [T5| . It reproduces the same 
reduced density matrix that satisfies a time-local master 
equation where the non-Markovian memory effect is fully 
taken into account. 

With the help of the exact reduced density matrix, the 
dynamics of an open quantum system could be well de- 
scribed. Meanwhile, there are several proposals to mea- 
sure the degree of the non-Markovianity of open quan- 
tum process (l9l . [20| . Very recently, the general non- 
Markovian dynamics of the environment on its surround- 
ing open quantum s yste m are explored within the ex- 
act master equation [21(. The question is how to probe 
the general non-Markovian dynamics. We thereby pro- 
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pose in this paper a promising approach to probe the 
time-dependent memory effect of a bath on a damped 
micro-cavity by coupling the cavity to a two-level atom 
dispersively. To probe the non-Markovianity of the dis- 
sipation of the single model EM field in a cavity, we let 
atoms of large detuning pass through the cavity. We 
found that the non-Markovianity of the bath is explicitly 
reflected by the atomic decoherence factor. In the week 
coupling region, the periodically reviving amplitude de- 
creases along with the cavity-bath coupling strength and 
decays to finally. On the the contrary, in the strong cou- 
pling region, the reviving amplitude increases with the 
coupling strength and almost does not decay in the ultra- 
strong coupling case, as a significant non-Markovian ef- 
fect [2l|. This atomic decoherence factor could be de- 
tected through the Ramsey interference in experiments. 

In the next section, we solve the Heisenberg equations 
of the unified quantum system plus bath model (Bose and 
Fermi) and obtain their formal solutions. In Sec. Ill, the 
derivation of the exact master equation of Bose system 
is presented. The exact master equation of Fermi case 
is addressed in Sec. IV. In Sec. V, we propose to probe 
the non-Markovian dynamics of a damped cavity with 
largely detuned two-level atoms. Finally, the summery 
of our main results is given in Sec. VI. Some detailed 
calculations are displayed in the Appendices. 



II. UNIFIED QUANTUM BATH MODEL AND 
FORMAL SOLUTION OF THE HEISENBERG 
EQUATIONS 

We consider an open quantum system S, which inter- 
acts with another large system B called bath. The com- 
bined system S + B is usually assumed to be closed, thus 
regarded as a Universe. The coupling of S to B will lead 
to the dissipation and dephasing of S. There are various 
types of bath, but the most commonly employed baths 
are modeled with non-interacting bosons and fermions. 
In this paper we consider the specific cases: a Bose sys- 
tem is surrounded by a Bose bath, or a Fermi system 
is immersed in a Fermi bath. Here, we first solve the 
Heisenberg equations for both the Bose and Fermi cases 
and obtain their formally exact solutions. 

The Universe Hamiltonian H = H s + Hf, + Hi nt is de- 
composed into three parts: the Hamiltonian of the sys- 
tem is taken to be a quadratic form 



bath is given by 



H K = 



M [ai,a 2 , 



j a N B 



(1) 



which describes N s linearly coupled bosons or fermions. 
aj(aj) is the annihilation (creation) operator of the ith 
mode of the system satisfying the commutation rela- 

= 8ti' (T corresponds to the boson and 



tion 



Oj, a. 



fermion, respectively) and M is a positive definite Her- 
mitian matrix. The Hamiltonian of the Bose or Fermi 



i=i 



(2) 



with the number of the uncoupled modes of the bath 
7Vf,(^> N s ) and annihilation (creation) operators bi 
(b,) which satisfy corresponding commutation relations 



bub), 



=5 a'. As proofed in the most usual envi- 



ronment coupled to the open system could be well ap- 
proximated as a collection of harmonic oscillators with 
linear quadratic couplings. Here, the interaction Hamil- 
tonian is taken as the form of 



N s N b 

H in t = Yl (^ a l b l + 7 lil b \ a i 

i=l 1=1 



(3) 



In the Heisenberg picture, the dynamics of the system 
is governed by the Heisenberg equations: 

d, (t) =-iJ2 M H a i (*)-*£ mbi (t) , (4) 

3 I 

b l {t) = -\w l b l {t)-iY,V>i{t). (5) 



For convenience, we introduce the (N s + Nj,) operator- 
valued vector 



c(t)= 



a, b 



[ai(t) , a 2 (t) o-nM , h(§ , b$ , • • ■ , b Nb (t)], 



and the (N s + Nb) x (N s + Nb) coefficient matrix 



M R 
i?t E 



(6) 



where 



R 



V21 V22 



VN S 1 VN s 2 



VlN b 
V2N b 



7 lN s N b 



and 

E = diag [wi, wa, ■ ■ ■ ,WN b ] ■ 
Then Eqs. <j4j) and ([5]) are re-expressed in a compact form 



-c(t) = -ms(t) 



(7) 



It follows from Eq. © that H is a time-independent 
Hermitian matrix. Consequently, the formal solution of 
Eq. ||7|| is given by 

c(t) = exp [-mt] c (0) ee U (t) c (0) , 
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where U (t) = cxp [— \Ht] is the time-evolution operator. 
Splitting the matrix IA (t) into four blocks 



U{t) 



[w(t)] NsXNs 

{P(t)] Nb ,N s 



N s xN b 



we obtain formal solution of Eq. ([7]) as 

a(t) =W(t)a(0)+T(t)b(0), 
b(t)=P(t)a{0) + Q{t)b(0). 



(8) 



(9) 
(10) 



The dynamics of total system is governed by these four 
time-dependent coefficient matrices W(t), T(t), Q(t), 
and P(t). 

Up to now, all the results are obtained by formal op- 
erations, since these coefficient matrices need to be de- 
termined by the differential equations. As shown in Ap- 
pendix B, there are some connections between these co- 
efficient matrices, which take a crucial role in derivation 
of the exact master equation. 



A. Differential equations of the coefficient matrices 

Substituting Eqs. © and ([TO]) into Eqs. (g]) and ©, 
we obtain the equations of the coefficient matrices 

W (t) = — i [MW (t) + RP (t)] , (11) 

f(t) = -i[MT (t) + RQ (t)] , (12) 

P(t) = -i[EP(t)+R*W(t)] , (13) 

Q (t) = -i [EQ (t) + i? f T (f)] , (14) 
with the initial conditions 

W (0) = I, T (0) = 0, P (0) =0,Q (0) = /. (15) 

Here, / is the identity matrix and is the null ma- 
trix. The differential equations of P (t) and Q (t) are 
integrated to yield 



P(t) = -i f dre-W*-^ I&W (t) dr, 
Jo 

-i f dre iEr R^T (r) + / 
Jo 



(16) 



(17) 



Then, we obtain the integrodifferential equations about 
W(t) and T (t): 

W (t) + \MW (t) + [ drG (t - t) W (t) = 0, (18) 
Jo 



t (t) + iMT (t) + / drG (t-r)T (r) 



iRe 



-iEt 



(19) 



Here the (N s x N s ) kernel matrix G (t) = Re-' lEt R) 
characterizes the non-Markovian memory structure of S. 
Defining the interaction spectral function 



we rewrite the element of the kernel matrix G (t) as 
Gij(t) = f dLuJ tJ {Lj)e- lut . 



Thus, the matrix G (t) is fully determined by the inter- 
action spectrum. 

On the other hand, the coefficient matrices W (t) and 
T (t) are not independent. By taking the Laplace trans- 
form of the integral differential equations (fT8| and (fT9]) , 
we get 



W[p] 
T[P] 



C(W)-- 
W[p}£ 



[p + iM - 
-iRe- iEt ) 



G{p)\ 



(20) 
(21) 



where £(••■) represents the Laplace transform. Conse- 
quently, after the inverse Laplace transform, the matrix 
T (t) is given by 



T (t) = -i / drW (t - r) Re 



-\Et 



(22) 



Thus the dynamics of S could be completely described 
by a single coefficient matrix W (t), It is well known that, 
under the Wigner-Weisskopff approximation, one can ob- 
tain the quantum Langevin equations of the operators 
of S by means of the approximate solution of Eqs. (fTO|) 
and (fTO|) together with the Heisenberg equations (0J and 
(|S]) Q. In this paper, it will be shown that the exact 
master equation of the reduced density matrix can also 
be obtained based on the formal solutions ([ MTOj) of the 
Heisenberg equations. And the Wigner-Weisskopff ap- 
proximation leads to the quantum Born-Markov master 
equation. 



III. BOSON CASE IN COHERENT-STATE 
REPRESENTATION 

In this section, we derive the exact master equation for 
N s coupled bosons in a Bose bath. In the Schrodinger 
picture, the total density matrix p(t) = U (t) p (0) J7 t (t) 
of S + B obeys the Liouville-von Neumann equation 
ihp (t) = [H, p (t)], where U (t) — cxp (—\Ht) is the time 
evolution operator of the total system. We assume that 
the total system is initially in the direct product initial 
state p (0) = p s (0) (8) pi, (0), with density matrices p s (0) 
and pb (0) of S and B, respectively. Through a lengthy 
calculation in Appendix C, the reduced density matrix 
of 5 is expressed in terms of the coherent state of the 
system 

p s (t) = / dp, (a, a) dp ff, f ) \a) (a\ 



K 



Ps (0) 



(23) 



with x = [xi, X2, ■ ■ ■ , xjv 8 ] (x = a, a', £, ^'). The propa- 
gator, which governs the dynamics of the reduced density 
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matrix, is defined as 
K 



with 



(at , a' , I ft ; t) = j dpi (z) (a, z\ U (t) 1 1) 

x(e' PbWC/H*)!"'^- (24) 



Here |z)(z = [zi,Z2, - ■ ■ ,ZN b }) is the coherent state of B. 
Different from the previous derivation fl3l - fl5| where the 
propagating function is obtained using the coherent-state 
path integral method and tracing over the environmen- 
tal degrees of freedom completely through the Feynman- 
Vernon influence functional, the propagator could also 
be evaluated in the coherent state representation by con- 
structing the explicit total wave function [f| 

C/ f (t) \a', z) = exp Is 1 (t) ■ a' + P (t) ■ z\ |0) , (25) 

as shown in Appendix C. It deserves to be noted that 
we have used the identities W (t) |0) = |0) and O (t) = 
W (t)OU(t). 



A. Propagating Function 

Generally speaking, the bath is initially in its thermal 
equilibrium state 



p b (o) = n 



i 



fi+i 



exp 



-0Eb 



(26) 



where /; = 1/ [exp (/3a;;) — 1] is the mean occupation 
number of the Ith bath mode at temperature T = 
1/(/cb/3)- In this case, the integral over the bath in the 
propagator (f24|) is carried out to give (please refer to Ap- 
pendix C for the detail), 



Ji(t)i+^jt(t)s' 



K(d\d',£,£\t)=A(t)exp 
where 

A^^dctlll + Vit)}- 1 ] 



■ (27) 



J 1 (t) = [I + V(t)]- 1 W(t), 
J 2 (t) = V[I + V(t)}-\ 

J 3 (t) = I-W* (t) [I + V{t)]- l W{t). 

This reproduces the propagating function obtained by 
the coherent state path-integral method in the previous 
works, e.g., Eq. (31) in [l3| . For convenience, we have 
introduced a new N s x N s Hermitian matrix V (t) = 
T(t)fT^(t). Utilizing the relationship in Eq. ([22} be- 
tween the matrices T (t) and W (t) , we have 



V 



(*)= I dn f dTzWirdGfo-TjWifo), (28) 
Jo Jo 



G(t)=Rfe- iEt R^. 



(29) 



Without any additional hypothesis, the exact propa- 
gating function of the reduced density matrix of S is ob- 
tained. The dynamics of S is governed by the single 
coefficient matrix W (t), which is determined by integral 
differential equation (fT5|) . And the influence of the bath 
on the dynamics of S is characterized by two memory- 
kernel matrices G (t) and G (t) [IM1- 



B. The Exact Non-Markovian Master equation for 
Bosons 



In the proceeding subsection, we have obtained the 
exact reduced density matrix of S as in Eq. (|23|) . Now we 
construct the master equation through its time derivative 



Ps 



dn (a, a') d^ ( ) \a) (a' 



t,?*;t) (^ Ps (0)||> 



(30) 



And it is found that the time differential of the propa- 
gating function takes the following form (please refer to 
Appendix D for the detail) 



^TKa' - Tr 



K-a 1 fr + iO + f) V a ,K 



-(V^k) fr-if2+f)a / +V^ff+2r)v a . J Pf, (31) 



with Hermitian matrices 
f (t) = 



(t) = v(t)-w (t) w- 1 (t) v(t)-v (w (t) w- 1 {tj) * , 

(32) 



r(t) = 

and 



W (t) W' 1 (t) + (W (t) IF- 1 (*)) T 



w (t) w- 1 (t) — (W (t) w- 1 (t) 



(33) 



(34) 



For coherent state defined in Eq. (|Al|) . there exist the 
following relations 

d \d) = a \d) , <Jt (d\ = (d\ eft , 

V Q \a) = a f \a) , V Q » (a| = (a\ a. 

With these mapping, we can construct the exact master 
equation of the reduced density matrix of the Bose system 
5, i.e. Eq. (32) in Ref. [3, 
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Ps (t) = -i [H s (t) , Ps (t)] + ]T [f y (t) + 2T tJ (t) 



OjPs 0) «! - ^Kp* (t) - ^p s ajc 



a\p s (t) dj - ^ajajps (t) - ip s (t) a ja \ 



(35) 



where i/ s = attla is the effective time-dependent Hamil- 
tonian of the system S. The diagonal elements of Cl (t) 
are the modified time-dependent frequencies of the differ- 
ent modes of S and the off-diagonals represent the new 
interaction strength between the modes of the system. 
Without Markov approximation, the dissipation of the 
system and the fluctuation of the bath could not be sep- 
arated. The original role of the bath is reflected by the 
time-dependent decay coefficients V (t) and f (t) [13| . 



C. From Wigner-Weisskopff Approximation to 
Markov Master Equation 

In this subsection, it will be shown that the Markov 
master equation can be obtained from the exact master 
equation by taking Wiger-Weisskopff approximation [9J, 
instead of making a direct Markov approximation [16| . 
Here the exact master equation is applied to the sim- 
plest dissipative system consisting of a single harmonic 
oscillator with frequency Slo an d a Bose environment. 
In this case, Cl, T(t), and T (t) are just time-dependent 
numbers instead of matrices, which are all determined 
by W{t) in Eqs. ([521153]) . Under the Wigner-Weisskopff 

I 



I 

approximation, the solution of Eq. (|18p is given by 

W(t) = cxp [-T t - i(fi + Aw)t] , (36) 

where 

T =TrJ(n ), (37) 

is the decay rate of the oscillator induced by the coupling 
to the vacuum and 



J (Oo) 

UJ — Qq 



dui, 



(38) 



is the small frequency shift, with the interaction spec- 
trum J (cj). It is easy to find that, in this case, the param- 
eters of the master equation become time-independent 



n = n a + Aoj, r = r , f = 2/ (n ) 7o , 



(39) 



where / (f2 ) is the mean occupation number of the os- 
cillator. As we know, T characterizes the dissipation of 
the system and T corresponds to the fluctuation of the 
bath. 

Then the Born-Markov master equation of a damped 
harmonic resonator is obtained as 



p s {t) = - i[ha)a, p s (t)} + [1 + /(fl )] T [2a Ps (t)a) - (atap s (t) + p s (t)a?a)] 
+ f{n )T [2aV,(t)a - (aa 1 p s {t) - p s (t)a<J)] . 



(40) 



It is known that, for a damped harmonic oscillator, the 
quantum Langevin equation of the number operator ob- 
tained from the Markov approximation is same as the 
one from the Wigner-Weisskopff approximation In 
this sense, these two approximations are equivalent. 



IV. FERMI CASE IN COHERENT STATE 
REPRESENTATION 



coherent state representation [23J, |24( reads 

Ps (t) = J dp (a, a') dp (£ C') |<3) (a'\ 

xK (a\a',C p t t) (|| Ps (0) |f) , (41) 

where the components of vectors a, a', and ^' are 
Grassmann variables, p s (0) is the initial state of S. And 
the initial state of the bath is still assumed to be the 
thermal state 



Pb 



(0)=IJ(l-/,)exp \-P$Eb 



(42) 



In the previous section, we obtained the exact master 
equation of the Bose system. Analogously, in the case of 
Fermi system, the reduced density matrix in the fermion 



where = 1/ [exp (/3u>i) + 1] is the mean occupation 
number of the lih Fermi mode with j3 = 1/ (kgT). Af- 
ter tracing over the degrees of freedom of the bath, we 
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find that the propagator is of the same form as the Bose 
case 1141 . 



K = A exp <3 f Ji|* + £'j\a! + aJ 2 a' + ^ J 3 £ 

but the matrices in K change into 

A(t) = det[[I-V(t)]- 1 ] , 



J 1 (t) = [I-V(t)]- 1 W(t), 
J 2 (t) = [I-V(t)]- 1 -I, 

J 3 (t)=W^ (t) [I-V^W^-I. 

After the same procedure as the Bose system, the exact 
master equation of the Fermi system is obtained as the 
same one given by Eq. (8) in [15| 



(i) = -i [&. (t) , Ps (t)] + £ [2T« (t) - fy (t) 



ajp s (i) aj - 2 a i a iPs (*) - ^ (*) a J a l 



(43) 



where if s (i), r(i) 
as the bosons' [hj 



and r (t) are defined in the same way 



V. PROBING NON-MARKOVIANIANITY OF 
AN OPEN QUANTUM SYSTEM 

In this section, we consider how to probe the non- 
Markovianity of a quantum dissipation process in a real- 
istic physical system. We understand that such an ideal 
probing scheme is usually based on the non-demolition 
measurement [25| . The interaction between the probing 
apparatus and the system to be detected commutes with 
the free Hamiltonian of the system, thus such kind of 
measurement does not change the energy of the system. 
But it will retain the information of the system on the 
probing apparatus. Such non-demolition interaction can 
be implemented in the cavity-QED as the dispersive in- 
teraction between the atom and cavity (26l. [27|. On the 
other hand, it is feasible to prepare and analyze a two- 
level Rydberg atom in a state corresponding to an arbi- 
trary point on the Bloch sphere in the quantum optics 
experiments. 

To realize the probing non-Markovianity in the cavity 
QED system, we consider an open quantum system: a 
single cavity mode coupled to its bath of many bosonic 
excitation modes resulting from the cavity leakage. Let 
an atom pass through the cavity, and then examine the 
quantum coherence of the atom. In this case, the atom 
could record the intrinsic information of the cavity field 
to accomplish the probing of the non-Markovianity of the 
cavity dynamics. This kind of approach was also used 
to p robe the quantum criticality of many body system 
|28| , where the sensitive change of the atom decoherence 
factor, which is characterized by the Loschmidt echo (29| . 
could reflect the quantum criticality of its surrounding 
environment. 




1 + ) 
I-) 



Figure 1. (Color online) Schematic diagram for probing of 
the non-Markovian dynamics of an open quantum system: a 
leaking cavity. The two-level atom passing through the cavity 
is largely detuned from the frequency of the cavity mode to 
approach the non-demolition measurement. 



In our case, the frequency of the atom uj a is drastically 
detuned from the the cavity resonance frequency ojq, i.e., 
A = luq — uj a 3> g a -fi where g a -f is the vacuum Rabi 
frequency characterizing the atom-cavity coupling. By 
making use of an adiabatic elimination procedure, we 
obtain the effective Hamiltonian 



fkuoa^a + hu a a z + hSa^c 



(44) 



for our probing scheme from the usual Jaynes-Cummings 
model |30|. Here, a(a') is the annihilation (creation) op- 
erator of the cavity, a z = |e) (e| — \g) (g\ is the Pauli 
matrix of the atom with the ground (excited) state of 
the atom \g) (|e)), and 8 = g\_fj A is the effective dis- 
persive coupling constant (3ll l-S2l | . Meanwhile, the cavity 
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Figure 2. (Color online) Decoherence factor with different 
cavity-bath coupling strength, (a) A = 0.002. (6) A = 0.01. 
(c) A = 0.1. 



Before entering the cavity, the atom is initialized in 
the superposition state (|e) + \g)) /y/2 and the cavity is 
initially in the coherent state \a). For simplicity, we as- 
sume that the bath is at zero temperature with initial 
density matrix p (0) = |0) (0|, where |0) is the vacuum 
state of the bath. It is well known that the bath of the 
cavity will decrease the coherence of the atom by disturb- 
ing the phase of the cavity field, but it does not change 
the population of the atom as the result of the disper- 
sive atom-cavity coupling. However, we can detect this 
decoherence effect by observing the Ramsey interference 
fringes of the out-coming atom. The exact density ma- 
trix of the atom and field is obtained by tracing over the 
degrees of freedom of the bath 



Pa-f 



= Tr 6 [t 



Mi 



(|^(0)) <^(0)|)®p 6 (0) 



iHtl 



(45) 



where H = H p + H b + fl int and \tp (0)) = (|e) + \g)) ® 
| a) /2. In order to describe the decoherence process of 
the atom, we introduce the decoherence factor |33| 



D{t) = \e 



I |2 

M Tr 



\Pa-f 



(46) 



where we have added a normalization factor cxp(— |a| ). 

If there were no bath present, the decoherence factor 
would read 



Do (t) 



1 



■ exp 



|2 I -Hit 



-1) 



(47) 



which is similar to the result in Ref. [32j. Thus the 
norm of the decoherence factor will decline to a very 
small value for |a| 2 3> 1 at the beginning and revive 
at St = nir,(n = 1,2,3,...) as depicted by the gray 
solid lines in Fig. 2. Since the cavity evolves along two- 
pronged path in the Hilbert space corresponding to dif- 
ferent atomic states and the two paths cross periodically. 

When the environment of the cavity is taken into ac- 
count, we obtain the decoherence factor from Eq. (|46|) 



D{t) = \e*v 



<w:{t)w al {t)+j^ a ,{t)-i)\c 



(48) 



where W a is determined by Eq. (|T8|) with M = ujq ± S (± 
corresponding to |e) and \g) states, respectively), and 



r t rt 

■h.aa'= I dT dT'W*,{T)W a {T r 



i) 







du)J(u))e 



;(t-t') 



Here, we choose the Ohmic spectral density with cut-off 
frequency O c : 



is coupled to a bosonic bath 

H b + fl^t = + {vi^h + H.c.) . 

i i 

Here, the atom has enough long coherence time and we 
neglect the decay of the atom during the strong probing 
process. 



J (us) = Awexp ( — — 



where A is a dimensionless constant characterizing cavity- 
bath coupling strength. 

Next we numerically calculate the norm of decoher- 
ence factor with or without Markov approximation with 
parameters: ujq = 1, 6 = 0.1, \a\ 2 = 5, and fl c = 10. 
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Figure 3. Norm of the decoherence factor without Markov ap- 
proximation, (a) If the cavity-bath coupling strength is weak, 
the recovering amplitude of the decoherence factor decreases 
along with A. (a) When the cavity-bath coupling strength is 
large enough, the recovering amplitude of the decoherence fac- 
tor increases along with A but its recovering period is changed 
by the bath. 



It is found that when the cavity-bath coupling is small 
(A = 0.002), the decoherence factors with or without 
Markov are nearly the same Figs. 2(a), but they di- 
verge from each other when the coupling strength be- 
comes large (A = 0.01) as in Figs. 2(b). And the Markov 
approximation loses its validity in strong-coupling regime 
(A = 0.1). From the insets of Figs. 2(a-c), we find that the 
Markov approximation also becomes invalid for a short- 
time dynamics (the norm of the decoherence factor under 
Markov approximation exceeds 0.5). 

When the cavity-bath coupling is weak, the decoher- 
ence factor without Markov approximation will still re- 
vive at St — nn,{n = 1,2,3,...), but the recovering 
amplitude decreases along with the cavity-bath coupling 
A and will decay to finally (Fig. 3 (a)), due to the 
dephasing of the cavity field induced by the bath. On 
the contrary, if the cavity-bath coupling becomes strong 
enough, the reviving magnitude will increase with the 
coupling strength A (Fig. 3(b)). Especially, when the 



Figure 4. Ramsey interference is used to detect the decoher- 
ence factor, (a) Real part of the decoherence factor in weak 
coupling region without the Markov approximation, (b) Real 
part of the decoherence factor in strong coupling region with- 
out the Markov approximation. 



coupling strength become to be ultra-strong (A > 0.1), 
the recovering amplitude almost does not decay, just 
like that the bath does not exist. This is because when 
A > wo/^c = 0.1 (for Ohmic bath), the cavity will stay 
in the system-bath coupling-induced dissipationless local- 
ized mode [2l[. As a result, the recovering amplitude al- 
most does not decay but the recovering period is shifted. 

Finally, we can utilize the Ramsey interference to de- 
tect the decoherence factor. After interacting with the 
cavity, the atom undergoes an additional resonant mi- 
crowave 7r/2 pulse performing the following transforma- 
tion 

|e>->^(|e> + |fl», |$>->^(-|e> + |s». 

And it is found that (please refer to Appendix E for the 
detailed calculation) 

Rc[D(t)] = ^[n g (t)-n e (t)}, (49) 
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where 

U a = e-^ 2 Tr f {{a\ e^' 2 p a . f e w ^l 2 \<r)} , (50) 

is the population of the atoms in the rotated state 
cxp(z0(jj,/2) |er) (tr = g, e) with rotation angle 9 = ir/2 
corresponding to the final tt/2 pulse. Thus we can mea- 
sure the real part of the decoherence factor through de- 
tecting the population difference of the out-coming atom. 
As shown in Figs. 4, the real part of the decoherence fac- 
tor can also reflect the non-Markovianity of the bath. 



where is the creation operator of the boson and \n) is 
the nth Fock state. It is found that the coherent state de- 
fined in Eq. ()A1[) is not normalized and different coherent 
states are generally not orthogonal 

(a | a') = (0| e Q * Q e Q ' at |0) = exp (aV) . (A2) 
All the coherent states form an over-complete sets 

J dn{a) | a) (a| = 1, (A3) 



VI. SUMMERY 

By constructing the reduced density matrix from the 
formal solution of the Heisenberg equations, we revis- 
ited the exact non-Markovian master equations for open 
quantum systems of Bose or Fermi type. The non- 
Markovianity can be reflected by the time-dependent de- 
cay coefficients such as T (t) and T (t) , with historical 
memory. To probe the non-Markovianity of the dissipa- 
tion of the single model EM field in a cavity, we let large 
detuning atoms pass through the cavity. It displayed 
that the non-Markovianity of the bath is explicitly re- 
flected by the atomic decoherence factor. In the week 
coupling regime, the periodically reviving amplitude de- 
creases along with the cavity-bath coupling strength A 
and decays to finally. However, in the strong coupling 
regime, the reviving amplitude increases with A and al- 
most does not decay in the ultra-strong coupling case. 
But the recovering period is shifted by the bath. We 
expect our results to be verified by experiments. 



with the measures 



dn (a) = e-H 2 — = e-H 2 -drdip. (A4) 

And the density matrix of the thermal equilibrium state 
in this coherent-state representation reads 

Pt = 1 - jr^: exp (-j3w a^a) (A5) 



dfx (a) 



1 



exp 



a 



/(wo) 



\a)(a\, (A6) 



where / (ujq) = 1/ [exp(/3cjo) — 1] is the mean occupation 
number, with temperature T = 1/(&b/3). 



2. Fermion coherent state 

The fermion coherent state is defined of a similar form 
as bosons 
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Appendix A: BOSON AND FERMION 
COHERENT STATES 

1. Boson coherent state 

For an arbitrary complex number a = rexp(i<p), the 
coherent state of a Bose mode with frequency ujq could 
be defined as 




The only difference lies in the fact that a is a generator 
of a Grassmann algebra instead of an ordinary complex 
number and is the creation operator for Fermi particles 
and they satisfy the anti-commutation relations 

{a,a'} = {a, a} = {a, a 1 } = 0. (A8) 

The overlap of two fermion coherent states is 

(a | a') = exp {a* a') , (A9) 

and the completeness relation reads 

J dfi{a) | a) (a\ = 1, (AlO) 

with 

dfj,(a) = da*dae~ a * a . (All) 



10 



Appendix B: Constrains of blocks of U (t) 

Due to hermiticity of matrix Ji, the time-evolution op- 
erator U (t) in Liouville space is a unitary matrix, i.e., 



W(t) T(t) 

Pit) Q(t) 



wHt) pt(t) 

Tt(t) Qt( t ) 



P 



which leads to 



PP f + QQ t = /, 



VFP f + TQ t = 0, 



fW f + QT+ = 0. 



(Bl) 



(B2) 



(B3) 



(B4) 



(B5) 



Except some special time t, the matrices W (t) and Q (t) 
are reversible. Then we have 

P= -QT^ (W^y 1 , (B6) 



P f = -W^TQl (B7) 



Here, we have used the fact the initial state of the total 
system is of the direct product form and the completeness 
of the coherent states {I-?)} and {|z"'}} of the bath. 

With the help of Eqs. ©, (TO]), (ESD, (El, and ([2~rJj) . 
the propagator is re-expressed in terms of the coefficient 
matrices 

K (eft, a',i,^;t) 
= J dp, (z) (w f a + P f z | t) (? I W^a' + 

x (TT-r^-) (T*5+QlZ\exA^I3Eb]\T&+QlA (C7) 

Using formulas (a\ cxp (Stfb) \a') = cxp [a*a' exp (5)] 
and 



<i/i (z) e 



cxp 



det [I - D] 



(C8) 



(for any N s x iV s Hermitian matrix D makes (/ — D) 
positive-definite), one goes to 



K(a\d',£,£ f ,t)=A(t)exp 
where 



+ ^ j\{t)a' 
+cflj 2 (t)c? + pj a g 



,(C9) 



Appendix C: CALCULATION OF THE 
PROPAGATING FUNCTION 



The reduced density matrix p s (t) of the system is ob- 
tained by tracing over the degrees of freedom of B in 
Pit) 



p s (t)= I dp{z*)(z\p{t)\z) (CI) 
' dpia, a)p s (<3, a'; t)\a)(a'\ , (C2) 



where \a) (|<5')) and \z) are coherent states of S and B, 
respectively. The element of the reduced density matrix 
is explicitly given by 



p s (d,a';t) = J dp(z)(d,z\p(t)\d',z) 



(C3) 



p s (0)p b (0)f, z"\tf\a', t) (C4) 

Jdp[C t)K(^,a', I f t ,*)(|jp s (0)||'} ) (C5) 



with 



K= j dp{z){d 1 z\Uit)\t) 

x(f p b io)uHt)\d\z) 



(C6) 



A = 



dot 



I-Qfil + f) x Q f , (CIO) 



J 1 =VF+T/(/+/)- 1 Qt / _Q / ( /+/ )-^gf Pj (cil 



J 2 =T/(/+/)- 1 Tt +Tf(I+f) 



/ - Q/ a+/) _1 Q t l Qfii+fV 1 t\ (Ci2) 



J3=P f 



(C13) 



and we have introduced a diagonal matrix / = 

diag[/i,/ 2l - •■ ,/jvJ- 

Then we will deal with these four terms one by one. 
First we make some pretreatment to obtain an expanding 
series. From Eqs. (|B3]|. (|B6|) . and (|B7)) . one finds 



/-Q/(/+/)- 1 Q t =Q[:rt(wt) V+Cj+^-^q^cu) 

So that 

i-Qfii+}V 1 Q ] 

oo 

=(gt)" 1 (/+/E(-i)"[r t (wt)" 1 r(/+/)]"g-. 1 (ci5) 

n=0 
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1. Ji (t), J 2 (t), and J 3 (t) 



In the second step, we carried out the integral over a of 
Eq. ([27| and used the identity 



According to Eqs. |B6|, (jCllj) . and (|C15j) . Ji (t) is 
explicitly expanded to 

oo 

r _ -i -in _ i 

= T^-T/^](-l) n p(WW t ) T{I+f) Tt(T^t) 

oo 

= W — T/Tt 5] (-i)i>w t )~ 1 T(/+/) rfl " (m^ 1 ')^ 1 

V [l+ (W r W rt ) _1 T(/+/)T t ^W 1 )" 1 



J 3 (t) + 4 (/ + V) Ji (t) = /. (C29) 
And the last step, the following formula is used 



dn {a') (a') n e Q '* Q = a r 



(C30) 



W 



W - V [WW ] +T(/+/)T t ] W 

(i + y) _1 i¥. 



(C16) 
(C17) 



The third step we have introduced a new N s x A^-matrix 
V (t) = T (t) /Tt (i). Similarly, one obtains 

J 3 =I- W 1 (1 + V)' 1 W. (C18) 
The calculation of J2 is a little more complicated 

j 2 =t|/ +fjr (-i) n frt(wt) -1 T(j+/)] n 1 

I n=0 J 

x/(J+/)- x Tt (C19) 

= y + t/ f (-i)"[rt(wt) _1 T(j+/)] n J 

xT* (W^'T/Tt (C20) 

(OO — 1\ 

xT t (W^ t )~ 1 T(/ + /-/)T t (C21) 
= V j 7 + £ (-ir[(^VKt)~ 1 T(/+/)Tt]"| 

00 

+^(-l)"[(lW t ) _1 T(/+/)Tt]"(MW t )^'rrt(C22) 
n=0 

= V [WW f +T(I + f) Tt] _1 (W+ + TTt) (C23) 



v(i + vo" 



(C24) 



2. A(t) 



The matrix A (i) is determined by the normalization 
condition, 

l = Tr[p s @] (C25) 
[<HcIHj(£, ?)K(a,a,g£, *)(l|p- (0) |f ) (C26) 

Adet[I+V]fd^, f )(e|p.(0) ^)exp(ft.|) (C27) 

Adet [I + V] /! £$/ |Y|k(0) le) • (C28) 



Since the initial density matrix is also normalized, thus 

A(t) = dot [(/ + vy 1 . 



Appendix D: TIME DIFFERENTIAL OF THE 
PROPAGATING FUNCTION 

The time differential of the propagating function is 
given by 



K 



A 



a f Ji I + ^ j\ a' + ^j 2 a'+ ft j 3 1* 



We define the differential operators 



and 



v„ = 



9 a 



(9 



A. (Dl) 



(D2) 



It is ready to find that 

£A = Jf 1 ( V«. - J 2 a' ) A', 



(D3) 



(D4) 



e n A= ( vl-tfJn) (4) 'a, 



ftAf^V^-atjVjA Jj-^V^-J^'jA. (D6) 



These relations lead to 



A = cJt 



Ka 



A 



Tr 



Ji ) J3J1 1 J2 



+a< 



JiJj 1 — J^yJA J3J1 1 



A 



V a .A 



V„,A 



•CM JsJf'V^.A 



4) M -1 ^ 



(D7) 



?tf Ka' - Tr 



A-a f (T + ifi + f ) V a .A 



- (v a JfVr- ift+f)a' + V Q , (f + 2r) V Q ,A, (D8) 
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with Hermitian matrices 



Appendix E: DECOHERENCE FACTOR 



r = v- ww~ x v - v [ww- 



ww- 1 + (ww- 1 ^ 



and 



= - 
2 



WW' 1 - ( ww- 



(D9) 



(D10) 



(Dll) 



The last step the following relations have been used 
~d 



dt 



(I+V) L {I+V)+(I+V) '(ww^p+v) (du) 

- (J + vy 1 v(i + vy 1 ] {i + v) 

+ + Vy 1 (WW' 1 ) (/ + V) (D13) 
-{I+Vy^V -(wW-^I+vt (D14) 



Through the approach in Appendix C, we can obtain 
the element of the reduced density in Eq (|45|) 



p a -f (a f ,a;a' f ,a) 



1 



A aa i exp 



_ +a*fJ2,<T<T'Ot'f + a*J 3 , acr >a 



(El) 



where in the case of zero temperature bath 

A aa . =1, J haa , = W a (t) , J 2iaa , = 0, J 3iaa , =Pl,P a . (E2) 

Here W a is determined by Eq. (TU) with M = UJ ±S (± 
corresponding to |e) and \g) states, respectively) and P a 
is given by Eq. (JTHJ). Following from Eq. (|4T)]) . we find 
that the population difference of the out-coming atom 
just gives the decoherence factor 



t . 



\J{\ j 3 jy=-{lW){^- 1 j-WW-\l+V)+V, (D15) 
and 

(4) 1 hJl l J2 



A 

A-* 



= — InA+Tr 
dt 



V [ww- 1 

-V 



ww-^v 



-Tr 



dt 

= -Tr T 



in (i+vy 



i-ii + vy 1 

Tr 



(D16) 



T + V(I+V) 1 (D17) 
(D18) 



u g (t) - n e (t) 



= e - |Q| "T r/ 



= D(t), 



-(e|e- ieCT */Va-/e i( W 2 |e) 



(E3) 
(E4) 



where = tt/2 and cr-y = i (\g) (e| — |e) (<?|). With the 
help of Eq. (|E1[) , we obtain the decoherence factor as 



D (tyRe{exp[$Y;(t)W a ,(t) + J 3 ^, (t) - l)\a\ 2 ]} . (E5) 
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